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We discuss the effects caused by a resonant cavity around a sample of a magnetic molecular crystal 
(such as Mni2-Ac), when a time dependent external magnetic field is applied parallel to the easy axis 
of the crystal. We show that the back action of the cavity field on the sample significantly increases 
the possibility of microwave emission. This radiation process can be supperradiance or a maser- 
like effect, depending on the strength of the dephasing. Our model provides further insight to the 
theoretical understanding of the bursts of electromagnetic radiation observed in recent experiments 
accompanying the resonant quantum tunneling of magnetization. The experimental findings up to 
now can all be explained as being a maser effect rather than superradiance. The results of our 
theory scale similarly to the experimental findings, i.e., with increasing sweep rate of the external 
magnetic field, the emission peaks are shifted towards higher field values. 

PACS numbers: 75.50.Xx, 42.50.Gy, 42.50.Fx 

I. INTRODUCTION 

Magnetic complex molecules have attracted a great deal of attention in recent years, because they have remarkable 
properties, related to their high magnetic anisotropy and large value of spinA* 2 ^: for Mn^Ac as well as for FesO, 
the quantum number of the total spin is S = 10. Accordingly, the eigenvalues of S z , (the spin component in the 
direction of the easy axis), can take 21 different values: rn = —10, . . . , 10. One of the most interesting properties of 
these molecules is that in a slowly changing external magnetic field the magnetization of the crystal consisting of such 
molecules exhibit series of steps at sufficiently low temperatures*^ The effect can be explained by assuming that the 
energy levels of the molecules become doubly degenerate at the corresponding values of the magnetic field and this 
resonance condition increases the possibility of the transition between the degenerate states with different values of 
m and m! . This kind of quantum tunneling between spin levels leads to a sudden change in the magnetic moment of 
the crystal and is therefore of fundamental importance as being a macroscopically observable quantum effect. 

In an important theoretical work Chudnovsky and GaraninS proposed that resonant magnetic tunneling could be 
accompanied by the emission of electromagnetic radiation, the possibility of superradiance from magnetic molecules 
has been considered further in Ref. [q and bursts of microwave pulses have actually been detected in recent 
experiments According to Refs. |"l6lll(| . the possible physical mechanism responsible for this phenomenon can 
be superradiance (SR)ii*i2iiii which is an interesting collective effect predicted first by Dicke in 1954, and has been 
experimentally observed in several physical systems since then^ However, when the radiation emitted by magnetic 
molecules was detected, the sample was placed in a container, which acted as a waveguide. This cavity changes the 
mode structure of the electromagnetic field surrounding the sample, which is known to have crucial effects on the dy- 
namics of the emitted radiation. Studies in SR with other physical systems like the ensemble of proton spinsi2iiiiSii& 
in the MHz, and with Rydberg atoms in the GHz domain^ show that the presence of a resonant cavity may en- 
hance the collectivity of the radiating individual dipoles, as first proposed by Bloembergen and PouncU^, and which 
seems to be necessary to obtain radiation in the case of molecular magnets, as welliA*!* 2 ^ Additionally, it has also 
been demonstrated that external resonators such as Fabry-Perot mirrors can enhance the relaxation of a crystal of 
molecular nanomagnets. 21 Inspired by these facts, we investigate in this paper the interplay between the radiation 
and the changes of the magnetization of a macroscopic sample of molecules Mni2-Ac inside a nearly resonant cavity. 
We also note that the interesting proposal of Ref. 22] to use these molecules for implementing a quantum algorithm 
gives another motivation to study their radiative properties. 

The present paper is organized as follows: In Sec.[H]we investigate the relevant magnetic level structure and describe 
a method that allows us to reduce the problem to a set of level pairs. The interaction of the molecules with the cavity 
field is considered in Sec. 11111 In Sec. llVl we discuss the approximate analytical consequences of our model and present 
numerical results as well. Finally we summarize and draw the conclusions (Sec. IVjl . 
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II. MAGNETIC LEVEL STRUCTURE 

Experiments including magnetization measurements neutron^ and EPR2£i2£i21 studies on crystals of Mn^Ac 
and Fe 8 suggest that the Hamiltonian responsible for the magnetic properties can be written as a sum of two terms: 

H S = H Q + H 1 . (1) 

Here Ho is diagonal in the eigenbasis of the (dimensionless) z component of the spin operator, S z : 

H = -DS 2 - FS A Z - flB S z , (2) 

where the last term describes the coupling to an external magnetic field applied in the z direction (easy axis): 
Bo = (0, 0, Bo) with ft = gfis- On the other hand, Hi consists of terms^*^ that do not commute with S z : 

H 1 = C(S% + St) + E(S\ + S 2 _)/2 + K(S+ + SJ)/2. (3) 

As most of the experiments where microwave radiation emitted by magnetic molecules was detected have been 
performed on Mni2-Ac, from now on, we shall consider this molecule as the representative example. In this case the 
values of the parameters in Ho are D/ks = 0.56K, and F/ks = 1-1 ■ 10~ 3 K. The coefficients in Hi do not have 
unanimously accepted values, but H\ can be considered as a small correction to Ho- However, as the transitions 
between levels with different m and m! are induced by terms that do not commute with S z , the importance of 
Hi is fundamental from this point of view. Tetragonal symmetry would allow only the quartic term, but there is 
experimental evidence^ showing the presence of weak quadratic and linear terms in Hi. We shall return to the 
determination of the coefficients in Hi at the end of this section. 

In a typical experimental situation the external field Bo slowly changes in time, and consequently so does Ho- 
Considering the total Hamiltonian as the generator of the time evolution (which means that relaxation effects are 
not taken into account), the corresponding time dependent Schrodinger equation governs the dynamics. However, it 
turns out that its solution is not feasible, because the time independent part of Hq forces a much faster evolution than 
the slow variation due to the change of the magnetic field. The fundamental frequencies w mm < ss (D/h)(m 12 — m 2 ) are 
in the range I0 10 — 10 11 s~ 1 being very fast compared with the time dependence of the magnetic field Bo that among 
ordinary circumstances^ cause a change on the scale 10s _1 . Therefore there is a significant variation in the molecular 
state as a consequence of the first term in @, while nothing happens due to —ftBoS z , which appears in the last term 
of Ho- On the other hand, we know that an appreciable change in the state occurs only if two energy levels become 
close to each other, as stipulated by a simple time dependent perturbation calculation, where the energy difference 
between the levels appears in the denominator of the transition probability. As the dominant term in the Hamiltonian 
is Ho, we can approximately find the points where two energy levels are close to each other by calculating the 
eigenvalues of Ho, which are obtained by the mere substitution S z — > m. Simple algebra shows that these eigenvalues 
become doubly degenerate with given m and m! at the following values of Bo'. 




A part of the level scheme of the total Hamiltonian H$ is shown in Fig. ^ as function of Bo. The special values of 
Bo given by Eq. |0J - where the levels of Ho cross - can be clearly identified in this figure. However, as it is known, 
the presence of Hi perturbs the eigenvalues leading to a splitting of the levels instead of crossing as shown by the 
inset. The resonance condition implies that appreciable changes in the population of the levels is expected around 
these avoided crossings (sometimes also called anticrossings) . Thus the system can be efficiently approximated by a 
set of level pairs, each of which is to be considered as an effective two-level system, similarly to the figure shown in 
the inset. 

Reducing the problem to a set of level pairs means technically that one applies degenerate perturbation theory 
around each avoided level crossing (determined by the values of Bo in Eq. Q) following a technique proposed 
independently by van Vleck2&, Des Cloiseaus2S and other authors as summarized in Ref. [3(J]. In the vicinity of 
a given avoided level crossing we find a unitary transformation that block diagonalizes Hs on the zero order two- 
dimensional eigensubspaces and generate an effective Hamiltonian H e that has the same eigenvalues as Hs'. 

H e = UH s tf (5) 

with 

U^P^Pay^P'aPa, (6) 

a 
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FIG. 1: The energy levels of Hamiltonian Q as a function of the external magnetic field Bo- The parameters are D/ttB — 0.56A', 
F/k B = 1.1 • W~ S K, K = 0.025/iBo, E/k B = -4.48 10" 3 K and C/k B = 1.36 10" 5 K. The inset zooms on the vicinity of 
Bo = 1.13 T, where - due to the presence of Hi - we can see a level splitting of Ao between the levels m = —10 and m = 8. 

where a labels the various eigensubspaces with quasi-degenerate eigenvalues, and P° are orthogonal projections on 
the degenerate eigensubspaces of Hq. For a given value of a, corresponding to the avoided crossing of levels m and 
m', we have P° = |m)(ra| + \m!)(m!\. As the perturbation is turned on, the {\m) , \m')} eigenstates of Hq evolve into 
{\4>m) , \<Pm'}} and we consider P = \4> m ){4> m \ + \4>m'){4>m' \ which projects onto the space arising from the zero order 
subspace of interest. By appropriate expansions^ii22^ of this unitary operator U , one can obtain a perturbation 
series, leading to useful analytical approximations of H e . Alternatively, the numerically exact projections P a can be 
found by the diagonalization of H$, and then Eqs. JSJ) and JBJ) provide the operator H e . The latter method is followed 
in this paper and we obtain in each two-dimensional subspace spanned by {\m) , \m')} the following matrix for the 
effective Hamiltonian: 

h m-f £o + w / 2 A °/ 2 ^ m 

H ^-{ A /2 eo-w/2 J , ( 7 ) 

with time dependent elements, and of course the values depend also on the pair {\m) , Here £o is the energy 

where the given crossing would occur, w is proportional to the time dependent external field in the z direction, while 
the offdiagonal clement Ao is the level splitting responsible for the effective coupling between the levels, w will be 
assumed to be linear in time with constant -Bo, yielding w(t) = —p,Bo(t — to)(rn — m') with to being the time instant 
when the crossing would occur. Note that this is a reasonable approximation even for a time scale much longer 
than the expected duration of the transition to be described (see e.g. Fig. 1. in Ref. 3). This linear approximation 
corresponds to the usual Landau-Zener-Stuckelberg (LZS) mode l 34 i 35 i 36 , by the aid of which we can calculate the 
probability of a given m — > m! transition: P m m' = 1 — exp(— 7tAq/2?iu;). Note that in this expression both w and Ao 
depend on the labels m,m'. For a given pair of levels and sweep rate Bq, P m , m ' is determined by the magnitude of 
the level splitting Ao, i.e., essentially by the parameters in Hi. If we assume that initially the system is in thermal 
equilibrium, we can consider a series of transitions at the values of Bq given by Eq. J2}. Calculating the expectation 
value of the operator S z (which is proportional to the magnetization) we obtain a staircase-like hysteresis loop that 
can be compared with the experimental curves^ 3 - at a given temperature and sweep rate. As our results depend on 
the coefficients in Hi, the minimization of the difference between the steps in the calculated hysteresis curve and 
the experimental plots gives the desired parameter values. We have found the best agreement for K — 0.025/ii?o, 
E/ks = —4.48 10~ 3 K, C/ks = 1-36 10~ 5 K, therefore these parameter values will be used in the following. 

The method summarized in this section, first of all, gives us information about the magnitude of the terms in the 
spin Hamiltonian, and describes how to obtain the level splittings Ao. As we shall see in the next section, the coupling 
of the molecular system to the cavity field at a given m — > m! transition can also be determined in this way. 
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III. INTERACTION WITH THE RESONANT CAVITY FIELD 



In this section we describe the interaction of an ensemble of magnetic molecules with a quasi-resonant cavity field. 
The time dependence of the level structure, considered in the previous section, will bring a certain level pair into 
resonance with the cavity field at a given value of the external magnetic field B . Usually there are avoided level 
crossings before the resonance, where for low lying m and ml the LZS transition probabilities P m .m' are small, leading 
to an almost complete inversion. We shall consider a dynamical equation for the density operator of the molecules 
and assume that a dipole moment is generated during a certain transition between the split levels m and ml of 
a given molecule which in turn serves as a source of microwave radiation influencing the transitions in all other 
molecules. To take into account this radiation mediated interaction, we have to add a new term to the effective 
two-level Hamiltonian describing the interaction with the magnetic dipole field of the cavity, B, which is an additional 
field beyond the stronger but almost static magnetic field Bo creating the inversion between the levels. The fast 
dynamics of the magnetic dipoles will be a forced oscillation generated by the interaction with the external field that 
is characterized by the interaction Hamiltonian Hj = jlBS e . Here S e = U SW — (S x , S y , S z ) denotes the spin operator 
one obtains after the unitary transformation JSJ) described in the previous section, by restricting it to the actual two- 
dimensional subspace we consider. The matrix elements of S e usually differ from those of S. As for the field strength 
B, a detailed model should take into account the mode structure of the cavity. However, the characteristic features 
of the dynamics due to the cavity field can be captured by a simpler model to be discussed here. We consider the 
microwave field as a single transverse (TM) mode being perpendicular to the z (easy) axis and having a frequency £1 
and equal amplitudes in the x and y directions: Hi — —jlB(S x + S y ) / \2. We note that the choice of the polarization 
does not have essential influence on the results presented here. Now the equations describing the dynamics of the 
two-level system (without relaxation) can be written as 

where g is the density operator of the effective two-level system and 

\ ' ' / mm' 

Here hu(t) = w(t) — 2fiBs' and A = Ao — 2jlBs with s' and s being the diagonal and offdiagonal elements resulting 
from the coupling operator (S x + S y )/v2. This leads to the following equations for the population differences and 
the coherences between the states: 

d i 

~^{Qmm Qm'm') — ~t Qmm' ^Qm'm) (10) 

d i 

~Jl8mm' = -iu{t)Q mm i + —A(Q mm - Qm'm')- (H) 

As is in the terahertz domain, we can separate a slowly varying amplitude of the time varying field and write 

B= (^B{t)e- mt +c.c\u k {zl (12) 

where Uk{z) is the corresponding mode function of the cavity, and B(t) <C f2 \B\. This makes straightforward a similar 

separation for the offdiagonal elements of the density matrix: g mm ' = R m m' e~ int where R m m' is again assumed to 
vary slowly compared with e~' int . Substituting into Eq. (|10|) we can neglect terms oscillating with frequency 2ft as 
they do not contribute essentially to the evolution of the state, i.e, this is the standard rotating wave approximation 
(RWA)iS. Introducing the notation Z mm i = g mm — g m 'm' for the inversion between levels m and ml , we obtain: 

d i 

~^ Z "irn' = —{(Rmm'^o e ' % ' _ R mm' ') ~ &{B* S* mm , R mm > — B S mm > R* mm ,)u\(z)} (13) 

d i 

-rRmm' — —i(u)(t) — n)R mm > + — (Aoe 1 — p,Bs mm >)Z mm > . (14) 
at n 

At a given time instant only one of the level pairs get into resonance with the cavity, therefore from now on we 
shall omit the indices m, m' . We shall come back to this point, and discuss the mechanism which selects the actual 
level pair. We can also average out the equations over a time period of a cycle of the oscillation that eliminates the 
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terms varying with frequency f2. A similar procedure can be performed in space over the wavelength A = 2ir/k, and 

we shall make use of j L u\(z) — 1/2. 

We also have to describe the effects caused by other degrees of freedom. These additional interactions - among 
which the strongest one is the spin phonon coupling, i.e, the oscillation of the atoms in the lattice - are not taken 
into account by the Hamiltonian . but can significantly influence the dynamics. The effects due to the reservoir of 
phonons (i) can be dissipative, by simply taking up the energy from the spin system and (ii) dephasing, by randomly 
disturbing the relative phases of the magnetic states. Dissipative terms lead to a decay of the diagonal elements while 
dephasing reduce the off-diagonal elements of g. The diagonal terms relax generally much slower than the offdiagonal 
ones, therefore we consider only this (so called transversal) relaxation. As usual, it will be taken into account by 
assuming a simple exponential decay with a time constant T2, the order of magnitude of which can be estimated 
between 10 -5 and 10 -7 s in the temperature range we are interested inS. With this term we have: 

±Z = -f n (B*s*R-BsR*), (15) 

^-R = -i(w(t) -Q)R- ^zfiBsZ - R/T 2 . (16) 
at 2n 

These equations are familiar from the theory of magnetic resonance, and as it is known, the effect of the field B is 
essential when the frequency originating from the slowly varying longitudinal field gets close to the cavity frequency: 
u(t) w fi. 

We also treat the mode amplitude of the cavity as a dynamical quantity. Following the usual semiclassical approach 
of radiation-matter interaction theory^ 7 -, the time varying field resulting from the magnetic molecules of the crystal 
will be described here as the field of a sample with time dependent magnetic dipole moment density M. in the x — y 
plane. The appropriate component of the transverse TL field originating from M. as a source, obeys the damped 
inhomogeneous wave equation 

AH-H/(c 2 T c )-H/c 2 = M/c 2 , (17) 

where T c is the cavity lifetime. Within the cavity we expand the field into modes, and in accordance with Eq. I|12|l . 
we also write: 

H=(±He- mt +c.c)ju k (z), M = (^-Me- int + c.c^f(z) (18) 

where f(z) is nonzero only within the sample, where it can be taken equal to Uk{z). If one substitutes into Eq. i|17|l . 
and makes an approximation exploiting that the amplitudes, H and M are slowly varying^ 7 , with respect to e~'°*, 
one obtains the equation 

dH „ M 1 

= iftn H. (19) 

dt ' 2 2T C V ' 

The filling factor 

77= / u k {z)f(z)dz/ [ u 2 k (z)dzfnl/L (20) 
Jc Jc 

arises when we project the resulting equation on the mode in question, by integrating over the volume of the cavity. 
Here l/L is the ratio of the lengths of the sample and the cavity, corresponding to the geometries reported in the 
experiments 7 ^. 

The corresponding component of the transverse magnetization of the sample is given by 



M = NofiTr 



g(S x + S y )/V2 =N il( S * ), (21) 



where Nq is the number density of the molecules participating in the transition m — > ml . We note that the static part 
of M. containing s' , does not give rise to radiation. Eq. I|21() connects the microscopic dynamical variables with the 
macroscopic ones. Recalling that Q mm < = Rmm<e~ int , we see that M = 2iV"o/2s*i? mm '. The slowly varying magnetic 
induction field acting on the molecules is given by B = [io(H + (3M), where /j,q is the vacuum permeability and (3 may 
differ from unity giving account of a local field correction resulting from the near field of the dipoles^. It is natural 
to measure the time variable in units of the characteristic time 

To= ™-^\ ■ (22) 
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As we shall see, the relation of To and the rate of relaxation characterizes the dynamics: in case of T0/T2 <C 1 the 
phase correlation of the individual emitters is conserved during the process and a superradiant pulse (or a sequence 
of pulses) can be emitted. On the other hand, To /T2 > 1 indicates that relaxation effects are too strong to allow SR 
to occur. 

To perform the calculations it is straightforward to introduce the dimensionless magnetic field strength and induction 
amplitudes: 

/ \ 1/2 / 1 \ 1/2 

h ={l5n) H > b ={j^m) B (23) 

The field intensity can be measured as the energy density averaged out over the time and space period. The di- 
mensionless field intensity / = |/i| 2 /2 gives the number of emitted photons of energy Ml per number of molecules 
participating in the given transition. Outside the sample h — b, while within the sample one has 

b = h + 2(3^-Re-^, (24) 

where ip is the phase of the off diagonal coupling constant: s — \s\e 1 ^ . The dynamical equation for the magnetic field 
in dimensionless form reads: 

^ = -l h + iRe-^, (25) 

where t = t/To and k — (To/T c ), is the damping coefficient of the cavity. These equations are to be solved together 
with Eqs. I|15llt)|) . which take the dimensionless form: 

^-Z = -i(b*Re-^ -WTe^), (26) 
dr 

-^-R = -iT (u(r)-Q)R-ibe**Z-'YR. (27) 

UT 

with 7 = T0/T2. Note that numerically (using SI units) To « | s | 1 x 10~ 8 s, where we substituted NqT] = 10 23 m -3 , 
corresponding to the values reported in the experiment^. Depending on the transition m — > m', usually the magnitude 
of the dimensionless matrix element \s\ is much less than unity, thus f Ons (obtained with \s\ = 1) is basically a lower 
bound for To ■ This value - at least at low temperatures - is less than the time scale of the relaxation, T2 , but clearly 
by orders of magnitude larger than the period of the microwave radiation thus our rotating wave approximation 
leading to Eqs. (|15lf 6fl is valid. Additionally, as is around 10 11 s _1 , Eq. I|24|) shows that it is a very good assumption 
to take h = b within the sample as well. 



IV. RESULTS AND DISCUSSION 



The generic scheme for the emission from the ensemble of magnetic molecules considered in this paper starts with 
inverted two-level systems that come into resonance with the cavity field at a certain value of the external magnetic 
field Bo- Besides resonance, an additional requirement for the transverse radiation to begin is that the wavelength 
corresponding to the transition frequency should be comparable or smaller than the size of the sample, otherwise the 
non-transverse near-field of the sources would dominate the emitted field at the location of the other molecules. This 
explains why in the experiments reported in Ref. [7J the emission is in the mm range wavelength. 

Now we shall analyze if the observed radiation can be considered as superradiance demanding 7 -C 1, or is it rather 
a maser effect, where the absence of phase relaxation is not crucial. Therefore we first assume that 7 = 0, and see 
that equations (|26|l and (|27() admit a simple constant of motion: 

Z 2 + 2\R\ 2 = Zl. (28) 

If lv(t) is changing sufficiently slowly, the condition of resonance lo(t) — Q c — is sustained during the dynamics of 
the emission. Then writing Z = Zq cos 9(t), and \R\ = (Zo/y/2) sin#(r), a simple equation yielding essential physical 
insight into the nature of the problem can be obtained. With the assumptions that R is real, = tt/2 and using 
Eqs. pfiTffT) . one has b = -§-6(t)/V2 = 6/V2, and from Eq. (gSJ we obtain 



0(r) + k9(t) /2 - Z Q sin 9(t) = 0, 



(29) 
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FIG. 2: Population difference Z (dotted line) and the dimensionless intensity (solid line) of the radiation emitted by the 
molecular system as a function of r for different cavity decay rates. The transversal relaxation rate 7 = and the dimensionless 
sweep rate corresponding to these plots is v — 0.2. Note that the energy per unit volume that corresponds to the field I — 1 
means that each active molecule in the sample emits a photon of energy hfl. 



which is the equation of a damped pendulum {9 measured from the inverted position), being often discussed in 
coherent atom-field interactions. The physically realistic initial condition for 9q is a small (6q ~ 0) value, as initially 
we expect the offdiagonal element of the density matrix to be small. This comes from the small initial polarization 
as a remnant of the coherence between the levels during the magnetic tunneling transition. Putting in such an initial 
condition takes into account the rapidly varying terms omitted in Eqs. I|14|) and (|13|l which would also lead to a small 
but nonzero R initially. Experimental results reported in Ref. 7\ suggest that the radiation appears after a crossing 
with small tunneling probability P mm >, thus there is a significant inversion present in the system, meaning Z ~ 1 in 
the beginning of the radiation process. 

If the cavity lifetime is short, k is large, the second term dominates over the first in Eq. (|29|l . Neglecting 8, the 
equation of the overdamped pendulum admits an analytical solution. For the emitted intensity one obtains 

\b\ 2 = 9 2 (r)/2 = -Vscch 2 [(r - r d )/r R ], (30) 

where tr = 2^' an( ^ Td = TR( m ^o/2). In this bad cavity limit T c can be estimated as L/c, the time needed for 
a photon to leave the cavity of length L. Assuming Zq = 1, the characteristic time of the emission in usual units 
is given by T7? = T$tr = Tq jj- — t)LAf ^~ p\s\' 2 ' ^ e see ^ me constant is inversely proportional to the 

number density of the molecules, ./Vo, while according to Eq. (|30|l . the intensity of the radiation is proportional to N§: 
these are the characteristic features of superradiancei^. In addition, there is a delay time T c i = t^Tq necessary for the 
appearance of the pulse described in Eq. (|30|l . 

In the case the cavity losses do not dominate the process, the energy of the field is fed back into the crystal and 
according to Eq. (|29[) this leads to a damped periodic process. Assuming a perfect cavity one has a kind of Rabi 
oscillations with a time dependent field: energy is exchanged periodically between the crystal and the field. 

These considerations based on the analytic solutions, however, become only valid approximately, as they do not take 
into account the factor u(t) — Q c in Eq. 127|) . As we assume a constant Bo, we can introduce a constant dimensionless 
external field sweep rate v via 

T (cj{t) - fi c ) = vt, (31) 

where the origin of the time axis is chosen so that r = corresponds to exact resonance: lo(0) = Q c . As an example, at 
the m = —10 — * to' = 8 transition with Bo — 30 mT/s and To = 10~ 6 s we have v = TojlBo{m — m')/h ss 0.1, leading 
to a dynamics significantly different from the analytical solution. Qualitatively, we expect that around t = (crossing 
point) the coupling begins to act, and creates a superposition of the levels. The coherence of the levels begin to increase 
accompanied by a finite transition probability: Z will differ substantially from Z{— 00). Then the oscillation of the 
pendulum, i.e., the radiation starts, but as the levels separate, their energy difference and therefore the oscillation 
frequency becomes larger. Taking relaxation into account, the amplitude of these oscillations diminishes, the molecules 
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FIG. 3: The effects of the transversal relaxation on the dynamics of the population difference Z (dotted line) and the intensity 
(solid line) of the radiation emitted by the molecular system. In comparison to Fig. [21 note the saturation effect in Z: for a 
given sweep rate (in this case v = 0.2) increasing values of 7 push the final population difference towards positive values. 

do not emit radiation any longer and simultaneously Z will reach a stationary value Z(oo), analogously to what is 
usually called quantum tunneling of magnetization, because a different Z means different value of the expectation 
value of S z . In this sense the inclusion of this time dependent detuning leads to a similar effect as discussed in the 
problem of tunneling, but the resonance condition is ensured by the inclusion of the time dependent cavity field, 
therefore the dynamics is more complicated than in LZS theory. 

Quantitatively, Fig. [21 shows the inversion and the intensity of the emitted radiation as a function of time for 
different strengths of the cavity decay. The effect of the resonance is clear, appreciable radiation and change in Z is 
seen after r = 0. As the cavity damping becomes stronger we have less oscillations in the emitted intensity and the 
process starts later. For small values of k the final inversion Z(oo) is determined by the sweep rate, but when cavity 
losses become significant, the energy of the molecular system is lost via the cavity field during the process leading 
to Z(oo) —1. A bad cavity (k 1) overdamps the pendulum and the system leaves the vicinity of the resonance 
before observable emission occurs. We note that the inversion Z(oo) is not necessarily closely related to the final 
magnetization of the sample, because following the photon emission, cascade (not purely two-level) transitions related 
to a given side of the two- well potential can also have a considerable probability. 

So far it has been assumed that the phase memory of the system is conserved, 7 = To /T 2 is small and accordingly 
the damping term — was neglected in Eq. (|27l) . This was the assumption that led to the coherent behavior of the 
molecules resulting in superradiant emission. However, in reality there are at least two reasons to consider nonzero 7. 
One of them is the spin-phonon coupling which is temperature dependent, thus can be reduced by cooling the sample. 
Additionally, if the size of the system is smaller than the wavelength the near field dipole-dipole coupling between 
the molecules becomes important and ca n be sh own to lead to an effective phase relaxation4i2£ Microscopic studies 
of the latter effect can be found in Refs. for a detailed recent work see Ref. 41]. At very low temperatures 

this effect can be even stronger than the homogeneous broadening mechanism caused by elastic collisions with the 
phonons. At temperatures around 2K, however, where the experiments observing the radiation have been performed, 
the dephasing is predominantly due to spin-phonon interactions instead of dipole-dipole coupling. 42-43 In the present 
work all these relaxation mechanisms are incorporated effectively by an appropriately chosen damping coefficient 7. 

The consequences of phase relaxation is shown in Fig. [21 where a moderate constant cavity decay (k — 0.2) is 
also taken into account. For weak dephasing, we have similar oscillations in Z and pulse structure as shown in 
Fig. [21 Increasing the value of 7 the coherent Rabi oscillations disappear. Additionally, the final inversion Z(oo) 
is a monotonically increasing function of 7, and this can be considered as a remarkable difference between the two 
decay mechanisms. Note that this saturation effect can be responsible for the additional steps in the hysteresis curve 
published in Ref. following the most pronounced one which is accompanied by microwave radiation: The nonzero 
population that remains on the upper level after the transition can lead to an observable change of the magnetization 
of the sample at a next avoided level crossing. 

In the case of strong dephasing, the time derivative of R can be neglected with respect of jR, and from Eq. I|27|l we 
obtain that R follows adiabatically the time dependence of b. Substituting back into Eq. I|26|) we obtain the following 
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vx (dimesionless B () ) 

FIG. 4: dM/dBo as a function of the dimensionless external magnetic field for different sweep rates. 7 = O.f , re = 1 and the 
dynamics is governed by Eqs. 1251261271 . 



rate equation description of the process: 



dr 7 2 + v 2 t 2 ' 

db K . ,„„. 
-r = -77 fo + — ■ (32) 

Here atomic coherence does not play any role, thus the process cannot be termed as superradiance, it is rather a 
maser, operating on the inverted magnetic levels. 

An important experimental result is that the position of the peaks in dM/dBo corresponding to the radiation 
process does depend on the external field sweep rate v. In our model this is related to the time spent by the sample 
around the resonance. For a slowly changing Bq, the dynamics is similar to the case of constant detuning, where an 
analytical solution is known, while increasing the value of v, an appropriate numerical solution of the dynamics is 
needed. In the case of superradiance, Fig. 0] shows dM/dBo as a function of vt, i.e., the dimensionless external field 
Bq. As we can see, for larger values of the sweep rate v, the height of the emission peak decreases and its position is 
shifted towards higher field values. The shift is in agreement with the experimental findings, while as a consequence 
of the coherent interaction, dM/dBo exhibits oscillations with sign changes. 

If we assume that the maser effect is responsible for the radiation and use the rate equations (|32|) to calculate the 
dynamics of the system, somewhat different results are obtained. As Fig. [S] shows, dM/dBo accompanied with maser 
radiation, scales similarly with v as in the case of weak dephasing: larger sweep rates correspond to peaks at higher 
Bo fields, thus taking the cavity effects into account, this scaling property is not characteristic for SR. However, the 
oscillations seen in the superradiant case are absent in Fig. |SJ 



We shall analyze now from the point of view of transversal relaxation if the observed radiation could be superra- 
diance. The reduction procedure summarized in Sec. ITU allows us to calculate the matrix element \s\ and thus the 
characteristic time (|22|) of the emission process for any transition. As the experimentally observed radiation peaks 
were around 1.4 T, we focus on this value of the external magnetic field. The level structure of the Hamiltonian Q 
provides the resonant transition frequency for a given transition m — > m! , as well as the population of the upper 
level according to the Boltzmann factor. In this way we can calculate To for any transition as a function of the 
temperature of the sample. As for SR to occur the dephasing rate 7 = X0/T2 must be small, so we should look for 
the transition with the minimal value of To. We found that below approx. 0.8 K the transition m = — 10 — * m! = 8 
provides the shortest To, while above this temperature the transition from m = — 6 to m! = 4 yields the minimal 
characteristic time. As Fig. [H] shows, for low temperatures, i.e., ground state tunneling, the minimal To is of the 
order of seconds. This is a consequence of the very small coupling coefficient s, corresponding to this transition. For 
higher temperatures when the transition m = — 6^m' = 4 provides the shortest characteristic time, To significantly 
decreases as a function of temperature. This is a consequence of the relation T cx l/(|s|yQVo) (see Eq. Q22fl). where 
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--• v=0.15 
v=0.20 
v=0.25 




1 2 3 

vi (dimesionless B ) 



FIG. 5: dM/dBo as a function of the dimensionless external magnetic field for different sweep rates. 7 = 1, k = 0.1 but in 
contrast with Fig. 2] the rate equations 1321 have been used to calculate the time evolution. 



-10->8 transition 
-6— >4 transition 




1.5 2.0 2.5 

Temperature (K) 



4.0 



FIG. 6: The minimal value (over the possible transitions) of To given by Eq. 1221 at Bo = 1.4 T, where microwave emission 
from the sample has been detected. The minimum at low temperatures corresponds to the transition m = — 10 —* m' — 8, 
while above approx. 0.8 K m = —6 — > ml = 4 provides the shortest To. 



iVo is temperature dependent. That is, above 0.8 K, the most favorable conditions for superradiance might be realized 
in the case of the transition in = —6 — * m! = 4, with a still strong temperature dependence of To. 

However, the energy emitted during a transition process is not necessarily the highest for the lowest To. In fact, 
the population of the m = — 6 level - which determines the maximum number of the active molecules - is not large 
enough to explain the magnitude of the emitted energy observed in a recent experiment. Ref. [j| reports on radiative 
bursts of duration of a few milliseconds, where (at 2 K) the total energy emitted by the sample was detected to be 
around 3 nJ. Using the parameters of the experiment^ we investigated all the possible transitions and found the best 
agreement with the experimental data for the transition m = — 8^m' = 6, giving a value of T) in the ms range and 
a total emitted energy to be around 1.5 nJ. (Note that for initial states below in = —8 the time scale of the process 
turns out to be too long, while for m > — 8 the number of active molecules is too small.) Thus our model predicts that 
the process having the most important role in producing the observed radiation is the transition m = —8 — ► in' = 6. 

As we have seen, the character of the emission depends on the ratio 7 = T0/T2, where T2 is decreasing with 
increasing temperature. According to Fig. El even for the shortest possible To and relatively weak dephasing, with T2 
around 10 -5 -10 -6 s (Ref. |42|). we have 7 > 1 at 2 K. The millisecond time scale obtained here and observed also 
in the experiments is clearly longer than the relaxation time T2. Thus - unless a yet unknown effect decreases the 
disturbance caused by phase relaxation - the process responsible for the experimentally observed 7 bursts of radiation 
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seems to be rather a maser effect than superradiance. 

The fact that no emission was seen for external fields lower than 1.4 T (where there can be a resonance as well) can 
be explained by the strength of the coupling to the transversal mode: \s\ is generally at least an order of magnitude 
larger for the possibly relevant transitions at 1.4 T than at previous resonances. However, in a good resonator - like 
the Fabry-Perot mirrors in Ref. |2l| - one may expect that several resonances have observable consequences, radiation 
and - as it has been detected - enhanced magnetic relaxation rates. Additionally, we note that at high sweep rates& 
the system passes not only a single cavity resonance during the emission process, and consecutive resonances can 
broaden the peaks in the dM/dB (B ) plots. 

V. SUMMARY AND CONCLUSIONS 

In this paper we developed a model for the interaction of a crystal of molecular magnets with the magnetic field of 
a surrounding cavity. The sample itself generates this transversal field B, while it also acts back on the molecules. 
The most important point of our treatment is that the cavity mode with fixed frequency £1 comes to resonance with a 
magnetic transition at a given value of the external longitudinal magnetic field. Around this resonance the interaction 
of the molecules with the mode significantly increases leading to an observable burst of electromagnetic radiation as 
well as a change in the magnetization of the sample. Our model can describe different mechanisms of this radiation, in 
fact, there is a continuous transition from superradiance to maser-like effects. The crucial parameter here is the ratio 
of two time scales, the characteristic time of the process and the dephasing time 7 = T0/T2. For small values of 7 the 
time evolution of the molecules is coherent allowing for the strong collective effect of superradiance in a cavity. In the 
case of strong dephasing, the sample still can emit electromagnetic radiation, but now the coherence of the molecules 
plays no role, the maser rate equations with time dependent detuning can describe the process. For moderate values 
of 7 we have a transition between the two processes. By calculating the intensity of the emitted radiation we have 
shown that with increasing the sweep rate of the external magnetic field, the emission peaks are shifted towards higher 
field values in accordance with the experimental results. This statement holds for both emission mechanisms, but the 
detailed functional dependencies arc different for SR and maser emission. 

Based on realistic approximations for To and T2, the process responsible for the experimentally observed bursts 
of electromagnetic radiation is most probably not superradiance, but rather a maser effect. The comparison of time 
resolved experiments on the emitted radiation with our theoretical results would provide the necessary information 
in order to settle this question. We expect that at very low temperatures, when spin-phonon relaxation is weaker, 
the collective features of the radiation may become dominant. While this is an interesting problem on its own, it is 
expected that the analysis of the radiated field can yield additional information on the process of quantum tunneling, 
as well as on the detailed properties of the interaction of these crystals and the field. 
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